Superconductivity of metallic boron in MgB 2 
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Boron in MgB2 forms layers of honeycomb lattices with magnesium as a space filler. Band struc- 
ture calculations indicate that Mg is substantially ionized, and the bands at the Fermi level derive 
mainly from B orbitals. Strong bonding with an ionic component and considerable metallic density 
of states yield a sizeable electron-phonon coupling. Using the rigid atomic sphere approximation 
and an analogy to Al, we estimate the coupling constant A to be of order 1. Together with high 
phonon frequencies, which we estimate via zone-center frozen phonon calculations to be between 
300 and 700 cm" 1 , this produces a high critical temperature, consistent with recent experiments. 
Thus MgB2 can be viewed as an analog of the long sought, but still hypothetical, superconducting- 
metallic hydrogen. 



Before the discovery of high-temperature supercon- 
ductors much effort was devoted to boosting the criti- 
cal temperature of conventional, BCS-Eliashberg super- 
conductors Q. An exotic and appealing idea going 
back to the early 60's, was that of metallic hydrogen 
The arguments were very simple: due to the light 
mass, the phonon frequencies in metallic hydrogen would 
be very high, of the order of several thousand K, and 
the prefactor in the BCS formula would be very large, 
so that even a moderate coupling constant would pro- 
vide a sizeable T c . This idea can be quantified as fol- 
lows: for monoatomic solids, the electron-phonon cou- 
pling (EPC) constant, A, which enters the BCS equation, 
can be written in the so-called McMillan-Hopfield form 
§, A = N(0) (I 2 ) /M(uj 2 ), where N(0) is the density 
of states (DOS) at the Fermi level per spin per atom, 
(/ 2 ) is the properly averaged electron-ion matrix ele- 
ment squared, M is the atomic mass and (w 2 ) is (again, 
properly averaged) the phonon frequency. The product 
M (w 2 ) does not depend on the mass, but on the force 
constants only Q, while 77 = N(Q) (/ 2 ) , also known as 
the Hopfield factor, is a purely electronic property. Cor- 
respondingly, light elements, everything else being the 
same, are beneficial for superconductivity. 

Lacking metallic hydrogen, attention was focused upon 
compounds with light elements: carbides, nitrides (ar- 
guably, the superconductivity in fullerenes was a discov- 
ery along this road). Indeed, many of them were 'high-T c 
superconductors' on the contemporary scale: 10-15 K. It 
was pointed out || that (l 2 ) is rather large in these ma- 
terials due to the relatively high ionicity (although not as 
high as in MgB2), but N(0) is rather small. This led to 
the suggestion of cubic MoN, which would have a larger 
N(0) than existing nitrides and carbides, as a hypothet- 
ical superconductor with T c > 30 K || . 

The recently discovered medium- T c superconductor 
MgB 2 Q with T c >39 K is clearly a continuation of the 
same idea. The main component, B, is even lighter than 
C and N. Furthermore, electronic structure calculations 
show that the compound is not only quite ionic with a 



reasonable DOS, but also has strong covalent B-B bond- 
ing (the bonding-antibonding splitting due to in-plane 
B-B hopping is about 6 eV) and thus exhibits strong 
electron-phonon interactions. Interestingly, unlike car- 
bides and nitrides, and similar to metallic hydrogen, elec- 
trons at the Fermi level (and below) arc predominantly B- 
like. Mg s-states are pushed up by the B p z orbitals and 
fully donate their electrons to the boron-derived conduc- 
tion bands. In the following we will describe the physics 
of such 'metallic' boron in detail, present an estimate of 
the EPC constant, and propose some routes for optimiz- 
ing T c in this kind of compound. 

MgB 2 occurs in the so-called AIB2 structure. Borons 
form a primitive honeycomb lattice, consisting of 
graphite- type sheets stacked with no displacement. The 
borons form hexagonal prisms with the base diameter 
of 3.5A nearly equal to the height. This creates large, 
nearly spherical pores for Mg. As in graphite, the intra- 
planar B-B bonds are much shorter than the distance be- 
tween the planes, and hence the B-B bonding is strongly 
anisotropic. However, the interplane bonds are only twice 
as long as the intraplane ones, as compared to the ratio 
of 2.4 in graphite, allowing for a significant interplane 
hopping. 

We have calculated the electronic structure of MgB 2 
using a general potential LAPW code 0. For the 
Rigid Atomic Spheres calculations we used the Stuttgart 
LMTO-TB code M . For the exchange-correlation poten- 
tial, the Generalized Gradient Approximation of Ref. || 
was employed. Despite the rather simple crystal struc- 
ture, very few electronic structure calculations for MgB 2 
have been reported (a model TB calculation of Burdett 
and Miller |n| and a recent full-potential LMTO study 
jn]]), and these have concentrated mainly on chemical 
bonding, paying hardly any attention to transport and 
electronic properties. The results of our LAPW calcula- 
tions are shown in Figs, [l] and 0. We note first that there 
is almost no valence charge inside the Mg MT sphere (less 
than 0.2 e). About half of the total valence charge resides 
inside the B spheres, and about the same amount in the 
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interstitials. This is partially due to the fact that the 
chosen LAPW setup employs rather small MT spheres 
for Mg. For the LMTO calculations we used an atomic 
sphere of nearly the size of the free Mg atom (up to 
3.13as), and obtained, as expected, a larger charge of 
2.8 electrons. However less than 25% of the charge has 
s-character. The remaining charge of p-, d-, and /- char- 
acter arises not from Mg electrons but rather from the 
tails of the B p-orbitals and contributions from the in- 
terstitials. In fact, one can say that Mg is fully ionized 
in this compound, however the electrons donated to the 
system are not localized on the anion, but rather are dis- 
tributed over the whole crystal. 

The resulting band structure can be easily understood 
in terms of the boron sublattice. We plotted in Fig. |l| 
the character of the bands (we show only the B p char- 
acter, since contributions from other orbitals near the 
Fermi level are very small) . We observe two B band sys- 
tems: two bands are derived from B p z states and four 
from B p x ,y All these bands are highly dispersive (light), 
the former being quite isotropic and the latter more two- 
dimensional. Both p z bands cross the Fermi level (in 
different parts of the Brillouin zone) , but only two bond- 
ing p XjV bands do so, and only near the T point. They 
form two small cylindrical Fermi surfaces around the T-A 
line (Fig. ||). However, due to their 2D character, they 
contribute more than 30% to the total N(Q). 

In contrast, the p z bands have 3D character, since the 
smaller intraplane distance compensates for a smaller 
(pp7r vs. ppc) hopping. In the nearest neighbor tight 
binding (TB) model their dispersion is £k = eo + 
2t ppa cos ck z ± tpp^y/Z + 2 cos aik+2 cos a 2 k+2cosa 3 k 
where ai,2,3 are the smallest in-plane lattice vectors. The 
on-site parameter eo can be found from the eigenvalue at 
the K point and is ~ 1.5 eV above the Fermi energy. We 
estimated t pprT and t ppT: from the LMTO calculations as 
~ 2.5 eV and ~1.5 eV, respectively. This model gives a 
very good description of the band structure near and be- 
low the Fermi, although the antibonding band acquires 
some additional dispersion by hybridizing with the Mg p 
band. The role of Mg in forming this band structure can 
be elucidated by removing the Mg atoms from the lattice 
entirely and repeating the calculations in this hypothet- 
ical structure. The in-plane dispersion of the both sets 
of bands at and below the Fermi level changes very little 
(for the pp~K bands there is hardly any change at all, while 
the ppa in-plane dispersion changes by ~ 10%). The k z 
dispersion of the p z bands is increased in MgB2 as com- 
pared with the hypothetical empty B2 lattice by about 
30%, and these bands shift down with respect to the p XtV 
bands by approximately 1 eV. This shift, as well as the 
additional dispersion, comes mainly from the hybridiza- 
tion with the empty Mg s band, which is correspond- 
ingly pushed further up, increasing the effective ionicity. 
Substantial k z dispersion of the p z bands produces the 
Fermi surface which is is approximately mirror-reflected 



with respect to a plane between the k z = and k z = tc/c 
planes, one pocket (electron-like) coming from the anti- 
bonding and the other (hole-like) from the bonding p z 
band. The two surfaces touch at one point on the K-H 
line and form a honeycomb tubular network, replicating 
in reciprocal space the boron lattice in real space. 

The resulting bands are fairly 3D: the average Fermi 
velocities are v x ,y = 4.90 ■ 10 7 cm/s and v z = 4.76 • 10 7 
cm/s. The plasma frequencies are u> px ,y = 7.1 eV and 
ui pz — 6.9 eV. Correspondingly, we predict fairly isotropic 
electrical resistivity, with the linear slope dp/dT « 
0.08A tr fifl cm/K. The total DOS at the Fermi level 
is N(0) = 0.36 states/spin-f.u. 

The above described band structure is typical for an 
sp— metal. What is not typical is that this particular 
sp— metal is held together by covalent bonding with a 
substantial ionic component, which inevitably leads to a 
strong electron-phonon interaction. In the following we 
present a semiquantitative estimate of the correspond- 
ing coupling constant, and argue that the fortunate com- 
bination of strong bonding, reasonable N(0), and high 
phonon frequency is responsible for the high transition 
temperature in this compound. 

In its most rigorous formulation the McMillan- Hopfield 
formula reads |lj] : 

A = n- 1 J2m a I j p) aj7S , (1) 

ij 

where the indices i , j run over all n ions in the crys- 
tal; a, j3 are Cartesian indices; (NIi a Ijp) is the electron- 
ion matrix element averaged over the Fermi surface; and 
&ia,j(3 is the standard force matrix. While this expres- 
sion is exact, a number of simplifications are needed to 
make it more practical. First, because of the large size of 
the Fermi surface, one can neglect the nondiagonal terms 
and reduce Eq. |l| to one unit cell Q : 

^E^V^&N" 1 ' ( 2 ) 

i i 

The quantity $u — d 2 E to t/dR 2 is a local quantity, 
which can be calculated from the total energy differ- 
ences of frozen zone-center phonons. The Hopfield fac- 
tor can be calculated in the Rigid Muffin Tin (or Rigid 
Atomic Sphere) Approximation |12] (RMTA), assum- 
ing that the change of the crystal potential due to an 
ion's displacement can be described by shifting the elec- 
tronic charge distribution rigidly inside the correspond- 
ing atomic sphere. Although expressions for 77 have been 
derived for arbitrary site symmetry Jill , these are com- 
plicated and we will use here simplified formulas JO] for- 
mally correct for cubic site symmetry. With this sim- 
plification, rf can be readily calculated from the LMTO 
potential parameters and partial DOS's. We show the 
results in Table 0. 
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In order to get an estimate of the phonon spectrum 
of the system, we calculated all zone-center modes, us- 
ing the full-potential LAPW method. There are four 
distinct modes juj: one silent mode, B\ g (two borons 
displaced along z in opposite directions), one doubly de- 
generate Raman mode, E 2g (in-plane displacements of 
borons), and two infrared-active modes, which do not in- 
volve changes of in-plane bonds: A 2u (B and Mg planes 
moving against each other) , and a doubly degenerate E\ u 
mode (B and Mg planes sliding along x, y). Their frequen- 
cies are, respectively, 690 cm -1 , 470 cm -1 , 390 cm -1 , and 
320 cm -1 , and their force constants are 390, 180, 70, and 
44 mRy/a^/atom. All calculations were carried out at 
the experimental lattice constant and c/a ratio 

Given the physics of the electronic structure described 
above, it seems likely that the lowest mode couples little 
with the electrons and that its softness is derived from its 
in-plane acoustic character. Thus, we excluded it from 
the calculations of A below. Correspondingly, the average 
inverse force matrix for boron is ($#) _1 = Mg 1 (u~ 2 ) = 

7.1a 2 3 /Ry, corresponding to {u~ 2 ) ~ 400 cm -1 . To- 
gether with the above value for 77 this gives A w 0.7 and 
the logarithmically averaged frequency of the same three 
modes is (690 • 470 2 • 390) 25 = (wi og ) w 500 cur 1 « 700 
K. We can now estimate the critical temperature accord- 
ing to the McMillan formula, T c = ^f^- exp[- 1.02(1 + 
A)/(A — /i* — /i*A)]. Using for the Coulomb pseudopo- 
tcntial \x* the commonly accepted value 0.1, we obtain 
T c w22 K. 

It is instructive to compare the calculation above with 
a typical sp superconductor, Al (Table |l|). RMTA no- 
toriously underestimates the electron-ion scattering in sp 
metals due to the weaker screening than in <i-metals. Also 
in Al, A is underestimated by a factor of 2. The screen- 
ing properties of MgB2 are similar to those of Al (they 
have similar Thomas- Fermi screening lengths), and the 
electronic properties are similar as well. It is therefore 
tempting to scale the calculated A for MgB2 by the same 
factor. The scaled results are also shown in Tab. || and 
the corresponding T c is approximately 70 K. Although 
this table looks impressively quantitative, these numbers 
should be regarded only as rough estimates. The RMTA 
in MgB2 is clearly a worse approximation than in Al. 
The boron site symmetry is far from cubic and, more- 
over, large differences in the atomic sphere radii lead 
to artificial potential jumps at the sphere boundary, a 



problem for which there is no remedy in RMTA. In other 
words, with regard to EPC, our calculations should be 
considered as a qualitative indication of a strong electron- 
phonon interaction. We can nevertheless be confident of 
our main qualitative conclusions. 

Our main conclusion is that MgB2, being essentially a 
metallic boron held together by covalent B-B and ionic 
B-Mg bonding, is electronically a typical sp metal with 
a typical DOS. Strong bonding induces strong electron- 
ion scattering and hence strong electron-phonon cou- 
pling. An additional benefit is the high frequency of the 
boron vibrations (while the force constants remain rea- 
sonably soft). Superconductivity is mainly due to boron. 
The light mass and correspondingly large zero-point vi- 
brations (> O.las) suggest a possibility of anharmonic 
and/or nonlinear EPC, and possible deviation of the iso- 
tope effect from 100%. Isovalent doping may be beneficial 
if it increases the density of states N(0). Lattice expan- 
sion due to Ca doping should lead to an overall increase 
of the density of states, and may provide the additional 
benefit of reduced p z -s-p z hopping. Another interesting 
dopant is Na, which should not only expand the lattice, 
but also decrease the Fermi level, exposing more of the 
Px,y bands, which may provide an additional contribution 
to A. 

Finally, let us outline the directions for further theo- 
retical investigation. First, EPC calculations beyond the 
RMTA {e.g., in the linear response formalism) are highly 
desirable and computationally feasible. Second, calcula- 
tions with full structure optimization for (hypothetical) 
CaB2 and BeB2, and virtual crystal calculations for Na 
doping should elucidate the effect of isovalent and hole 
doping, giving some hints toward further optimizing T c . 
Work along these lines is currently in progress. 

We are thankful to J.E. Pask, OS. Hellberg and D.J. 
Singh for critical reading of the manuscript. This re- 
search was supported in part by ONR and by DOE under 
Contract No. W-7405-82. 

TABLE I. Partial LMTO DOS, N(0), in eV _1 /spin and 
partial Hopfield factors, 77, in mRy/a%, for Mg and B (per 
atom) . 
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TABLE II. Comparison of the electron-phonon coupling in Al and in MgB2. The entries labeled by a question mark are 
obtained by scaling the RMTA A by the ratio of \ exp /^rmt in Al. qrF is the Thomas- Fermi screening parameter in a^ 1 . T c was 
calculated by the McMillan formula with /j* = 0.1. LAPW iV(0) is given in eV _1 /spin, 77 in mRy/o| and (^Muj 2 ^ in Ry/a%. 
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FIG. 1. Bandstructure of MgB2 with the B p-character. 
The radii of the red (black) circles are proportional to the B 
Pz (B p x , y ) character. 
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FIG. 2. Total density of states (DOS) and partial DOS for 
the MgB2 compound. The small Mg DOS is partially due to 
the small rMT of 1.8 og used. 




FIG. 3. The Fermi surface of MgB2. Green and blue cylin- 
ders (hole-like) come from the bonding p x , y bands, the blue 
tubular network (hole-like) from the bonding p z bands, and 
the red (electron-like) tubular network from the antibonding 
p z band. The last two surfaces touch at the K-point. 
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